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ABSTRACT 



New H-band Galaxy Number Counts: A Large Local Hole 
in the Galaxy Distribution? 
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0\ ' We examine ii-band number counts determined using new photometry over two 

fields with a combined solid angle of 0.30 deg 2 to H w 19, as well as bright data 
(H ^ 14) from the 2 Micron All Sky Survey (2MASS). First, we examine the bright 
number counts from 2MASS extracted for the «4000 deg 2 APM survey area situated 
around the southern galactic pole. We find a deficiency of «25 per cent at H = 13 with 
^[ ' respect to homogeneous predictions, in line with previous results in the B-band and 

00 , _fsT s -band. In addition we examine the bright counts extracted for |6| >20° (covering 

Q\ ■ «27 000 deg 2 ); we find a relatively constant deficit in the counts of «15-20 per cent 

(^*> ' to H — 14. We investigate various possible causes for these results; namely, errors 

. in the model normalisation, unexpected luminosity evolution (at low and high red- 

shifts), errors in the photometry, incompleteness and large-scale structure. In order 
(-h \ to address the issue of the model normalisation, we examine the number counts de- 

O termined for the new faint photometry presented in this work and also for faint data 

(H < 20) covering 0.39 deg 2 from the Las Campanas Infra Red Survey (LCIRS). In 
each case a zeropoint is chosen to match that of the 2MASS photometry at bright 
magnitudes using several hundred matched point sources in each case. We find a large 
offset between 2MASS and the LCIRS data of 0.28±0.01 magnitudes. Applying a 
consistent zeropoint, the faint data, covering a combined solid angle of 0.69 deg 2 , is 
in good agreement with the homogeneous prediction used previously, with a best fit 
normalisation a factor of 1.095^0 034 higher. We examine possible effects arising from 



. unexpected galaxy evolution and photometric errors and find no evidence for a sig- 

nificant contribution from either. However, incompleteness in the 2MASS catalogue 
(< 10 per cent) and in the faint data (likely to be at the few per cent level) may have 
a significant contribution. Addressing the contribution from large-scale structure, we 
estimate the cosmic variance in the bright counts over the APM survey area and for 
|b| >20° expected in a ACDM cosmology using 27 mock 2MASS catalogues constructed 
from the ACDM Hubble Volume simulation. Accounting for the model normalisation 
uncertainty and taking an upper limit for the effect arising from incompleteness, the 
APM survey area bright counts are in line with a rare fluctuation in the local galaxy 
distribution of m 2.5o\ However, the \b\ >20° counts represent a 4. Oct fluctuation, and 
imply a local hole which extends over the entire local galaxy distribution and is at 
odds with ACDM. The increase in faint near infrared data from the UK Infrared Deep 
Sky Survey (UKIDSS) should help to resolve this issue. 

Key words: galaxies: photometry - cosmology: observations - large-scale structure 
of the Universe - infrared: galaxies 



1 INTRODUCTION 

A recurring problem arising from the study of bright galaxy 
number counts has been the measured deficiency of galax- 
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ies around the so uthern g a lactic pole. This was first exam- 
ined in detail by [Shank 4 dl990l) a nd sub sequently by the 
APM galaxy survey Jlvladdox etaTlll990all . which observed 
a large deficit in the number counts («50 per cent at B =16, 
~30 per cent at B —17) over a ~4000 deg 2 solid angle. If 
this anomaly was due solely to features in the galaxy distri- 
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bution, this would be at odds with recent measurements 
of the variance of l ocal galaxy densi t y fluctuations ( e.g. 
lHawkins et aDl2003l: iFrith et alJhoOSbl : Icole et alJkoOSI) or 
the expected linear growth of density inhomogeneities at 
larg e scales. 

iMaddox et alJ ljl990bl) examined possible causes of this 
deficiency . From redshift surve y results over the APM sur- 
vey area ijLovedav et alJll992l) . it was argued that a weak 
local under-density contributed to the obse rved deficiency 
at the < 10 per cent level at B ~17. Instead. IMaddox et alJ 
( 1990b) suggested that strong low redshift galaxy evolution 
was the dominant contribution. This phenomenon has also 
been suggested as a possible explanation fo r large deficien - 
cies in the Sloan Digital Sky Survey fSDSS: lLovedavll2004|) . 
although models without such strong low redshift evolution 
provide predictions co nsistent with observed number red- 
shift distributions (e.g.|B roadhurst et al. 1988; Collcs set alJ 
Il990l: lHawkins et aLlfeoofl l. In contrast. IShanksl l|l990l) ~ 
gued that evolution could not account for the observed slope 
and that large-scale structure was the principal cause of the 
deficiency in the counts. 

However, another possible contribution to the low 
counts might be errors in the APM photometry. Comparing 
the photograph i c APM photometry with B-band CCD data, 
iMetcalfe et, alJ l)l995l) detected a small residual scale error 
in the APM survey zeropoints for B> 17. Correcting for this 
offset, the counts were now in good agreement with homoge- 
neous predictions at faint magnitudes (B>17.5); however, 
the problemat ic deficiency at brighter magnitudes remained. 
More recently. IBusswell et al.l (|200j) used B-band CCD data 
over «337 deg 2 within the APM survey area to provide the 
most accurate comparison to date with a sample of the APM 
survey photometry. The photometric zeropoint of this CCD 
data was in excellent agr eement with the Millennium Galaxy 
Catalogue (Driver 2003) and the Sloan Digital Sky Survey 
Early Data Release (lYasuda et al.ll200ll) . However, a com- 
parison with the APM photometry suggested a large off- 
set of 0.31 magnitudes for B < 17.35. Applying this to the 
APM surve y counts, a defi c iency of «25 per cent remained 
at B =16; IBusswell et alJ i2004) determined that such a 
deficiency in the local galaxy distribution would still be at 
odds with a ACDM form to the galaxy correlation function 
and power spectrum at large scales. 

In order to examine this issue independently, bright 
numbe r counts have also been exam ined in the near in- 
frared llFrith et al.ll2003l. 12004 l2005aT) . These wavelengths 
are particularly useful for such analysis as the number 
count predictions are fairly insensitive to the evolution- 
ary model or the assumed cosmology at bright magnitudes 
(see Fig. 0; current observations are in remarkable agree- 
ment with predictions in the A-band to K ~23 for ex- 
ample jMcCracken et all |2000T) In particular, IFrith et al] 
(2005a) examined A s -band number counts select ed from the 
2 Micron All Sky Survey (2MASS: Ijarretdl2004f) . First, the 
counts over the APM survey area were determined; a simi- 
lar deficiency was observed to the APM survey counts ( with 
the zeropoint offset determined bv lBusswell et al] j2004t) ap- 
plied), with a ~25 per cent deficit at K s =12 co mpared to 
the no evolution model of IMetcalfe et all (1200 ll) . Using a 
ACDM form for the angular correlation function at large 
scales and assuming the observed counts were solely due 
to features in the local galaxy distribution, the observed 



counts represented a 5a fluctuation. However, this result 
was complicated by the fact that the 2MASS AVband num- 
ber counts for almost the entire survey (|6| >20°, covering 
f»27000 deg 2 ) were also low, with a constant deficiency of 
~20 per cent between K 3 =10 and K s =13.5. 

Did this surprising resu lt perhaps indicate that the K 3 - 
band IMetcalfe et alJ i200 II) model normalisation was too 
high? Or, as suggested previously, could low redshift lumi- 
nosity evolution significantly affect the br i ght cou nts? These 
issues were also addressed bv IFrith et all ll2005aD : First, the 
IMetcalfe et alJ (1200 ill model was compared with faint K- 
band data collated from the literature. Fitting in the mag- 
nitude range 14< K <18 it was found that the best fit 
model normalisation was slightly too high, although not sig- 
nificantly (this magnitude range was used so as to avoid 
fluctuations in the counts arising from large-scale structure 
at bright magnitudes and significant effects from galaxy 
evolution at the faint end). Accounting for the normalisa- 
tion uncertainty (of ±6 per cent) the observed deficiency 
in the A s -band counts over the APM survey area still rep- 
resented a ps 3cr fluctuation. Second, the issue of low red- 
shift luminosity evolution was also addressed: 2MASS galax- 
ies below K a = 13.5 were matched with the Northern and 
Southern areas of t he 2dF Galaxy Redshift Survey (2dGRS; 
IColless et aT] 1200 ll) . The resulting n(z), covering > 1000 
deg 2 in total, was con sistent with the no evolution model of 
IMetcalfe et alJ i200ll) . In addition, these A s -band redshift 
distributions were used to form predictions for the num- 
ber counts over the Northern and Southern 2dFGRS areas 
respectively. This was done by multiplying the luminosity 
function parame ter 4>* (which governs the model normalisa- 
tion) used in the lMetcalfe et all i200ll) model by the relative 
density observed in the A s -band n(z) as a function of red- 
shift. These 'variable (f>* models' were then compared with 
2MASS counts extracted for the 2dFGRS areas in order to 
determine whether the observed counts were consistent with 
being due solely to features in the local galaxy distribution; 
the variable <j>* models were in good agreement with the 
number counts, indicating that low redshift luminosity evo- 
lution is unlikely to have a significant impact on the observed 
deficiency in the counts, in the A s -band at least. 

In this paper we aim to address the issue of low, bright 
number counts in the near infrared ff-band. In particular 
we wish to address a drawback to the A s -band analysis of 
IFrith et alJ J2005al) - the issue of the number count model 
normalisation; while the A s -band model used was compared 
with faint data and was found to be in good agreement, 
the level to which systematic effects, arising perhaps via 
zeropoint offsets between the bright and faint data or cosmic 
variance in the faint data, might affect the conclusions were 
uncertain. We address this issue in the if-band using new 
faint data covering 0.3 deg 2 to H = 18, calibrated to match 
the 2MASS zeropoint. In section 2, we first verify that the 
_ff-band counts provide number counts over the APM survey 
area which are consistent with the previous results in the B 
and AVbands llBusswell et al.l2004l:lFrith et al]l2005al) . and 
that the form of the counts is not significantly affected by 
low redshift luminosity evolution through comparisons with 
the variable <f>* models described above. In section 3, we 
provide details of the data reduction of the new faint H-band 
photometry. The associated counts are presented in section 
4. In section 5 we discuss possible systematics affecting the 
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Figure 1. H-band galaxy number counts collated from the liter- 
ature. The dashed and solid lines indicate the no evolution and 
pure luminosity evolution predictions respectively, described in 
section 2. We also show bright H-band counts extracted from the 
2MASS extended source catalogue for \b\ >20°. For each dataset, 
we indicate the associated observed solid angle in square brackets. 



bright number counts including the model normalisation and 
incompleteness. The conclusions follow in sections 6. 



2 BRIGHT H-BAND COUNTS FROM 2MASS 

We wish to examine the form of bright number counts in 
the H-band in order to verify that the counts over the APM 
survey area («4000 deg 2 around the southern galactic pole) 
are comparable to those measured pre viously in the optical 
H-band and near infrared JG-band iBusswell et alJ 120041: 
iFrith et alJl2005al) . The near infrared has the advantage of 
being sensitive to the underlying stellar mass and is much 
less affected by recent star formation history than optical 
wavelengths. For this reason, number count predictions in 
the near infrared are insensitive to the evolutionary model 
at bright magnitudes. In Fig. Q we show faint H-band data 
collated from the literature along with bright counts ex- 
tracted from 2MASS over «27000 deg 2 . The 2MASS magni- 
tudes are determined via the 2MASS H-band extrapolated 
magnitude; this form of magnitude estimator has previously 
been shown to b e an excellent est imate of the total flux 
in the 7f s -band iFrith et alJ l2005bO t hrough comparison 
with the total magnitude esti mates of | Jones et al.l (I2004D 
and the A"-band photometry of lLovedaTl200ol) . Through- 
out this paper we use 2MASS H-band counts determined via 
this magnitude estimator. We also show two models in Fig.0 
corresponding to homogeneous predictions assuming no evo- 
lution and pure luminosity evolution models. These are con- 
structed from the H-band lum inosity function parameters 
listed in iMetcalfe et alJ (120051) and the K + E- correct ions 



Figure 2. Number redshift histograms for 11501 and 13 687 
H <14 2MASS galaxies matched with the 446 deg 2 Northern 
(left hand) and 647 deg 2 Southern (right hand panels) 2dFGRS 
declination strips respectively. In each case the solid lines indicate 
the pure luminosity evolution prediction for a homogeneous dis- 
tribution described in section 2 normalised by the respective solid 
angles. We also indicate the relative density in the lower panels, 
dividing the observed n(z) by the homogeneous prediction. 



of iBruzual fc CharlolJ lll993l) . At bright magnitudes the two 
are indistinguishable; only at H> 18 do the model predic- 
tions begin to separate. The faint data is in good agreement 
with both the no evolution and pure luminosity evolution 
predictions to H «26. 

Before examining the H-band counts over the APM sur- 
vey area, we first verify that the bright counts are consis- 
tent with relatively insignificant levels of low re dshift lumi- 
nosity evolution in the manner carried out by IFrith et alJ 
j2005al) for the H s -band counts. In the upper panels of 
Fig.|5|we show H-band n(z) to the 2MASS limiting magni- 
tude of H = 14, determined through matched 2MASS and 
2dFGRS galaxies over the 2dFGRS Northern (left hand) 
and Southern (rig ht hand panels) declination strips (see 
IFrith et alJ i2005al) for further details of the matching tech- 
nique). The solid lines indicate the expected homogeneous 
distribution c onstructed from the p ure luminosity evolution 
predictions of lMetcalfe et all i2005h (there is no discernible 
difference between this and the no evolution prediction). In 
the lower panels we divide through by this prediction; these 
panels show the relative density as a function of redshift. 
The observed n(z) are consistent with the expected trends, 
with relatively homogeneous distributions beyond z = 0.1 
(1 per cent and 8 per cent over-dense in the North and 
South respectively for 0.1 ^ z ^ 0.2). For this reason, Fig. [5] 
suggests that the level of luminosity evolution is relatively 
insignificant at low redshifts in the H-band; strong luminos- 
ity evolution produces an extended tail in the predicted n(z) 
which is not observed in the data. 




Figure 3. //-band 2MASS galaxy number counts extracted from 
the Northern (left hand) and Southern (right hand) 2dFGRS dec- 
lination strips. The solid line indicates the homogeneous pure 
luminosity evolution prediction described in section 2 (this and 
the no evolution prediction are indistinguishable at these mag- 
nitudes). The dashed and dot-dashed lines indicate the variable 
(ft* models for the Northern and Southern 2dFGRS strips respec- 
tively; these indicate the expected number counts given the ob- 
served n(z) (Fig.|2j. In the lower panels we divide through by the 
homogeneous prediction. In each case the errorbars indicate the 
Poisson uncertainty in each bin. 



Figure 4. _ff-band 2MASS galaxy number counts extracted for 
the APM survey area (Ri4 000 deg 2 ) and for |6| >20° («27000 
deg 2 ) , shown in the left and right hand panels respectively. As in 
Fig|5] we show the homogeneous pure luminosity prediction (solid 
line), and the Northern (dashed) and Southern (dot-dashed) vari- 
able <j>* models, indicating the expected number counts for the 
redshift distributions shown in Fig. [3] As before, in the lower 
panels we divide through by the homogeneous prediction. In each 
case the errorbars indicate the Poisson uncertainty in each bin. 



As a further check against strong low redshift luminos- 
ity evolution, we can use the observed n(z) to predict the ex- 
pected //-band number counts over the 2dFGRS declination 
strips . This technique is described in detail in iFrith et alJ 
J2003tl2005ar> . To recap, we use the observed density (Fig.|21 
lower panels), to vary the luminosity fu nction normalisation 
(<f)*) used in the Metcalf e et al.l (l2005t) model as a function 
of redshift (for z ^ 0.2). We show these 'variable <f>* models' 
along with the 2MASS //-band counts extracted for the 2dF- 
GRS strips in Fig. [3] In each case, the upper panels indicate 
the number count on a logarithmic scale; in the lower panels 
we divide through by the homogeneous prediction. In both 
the Northern and Southern 2dFGRS areas, the counts are 
in good agreement with the expected trend, defined by the 
corresponding variable (j>* model. This indicates that real 
features in the local galaxy distribution are the dominant 
factor in the form of the observed if-band number counts, 
and that strong low redshift luminosity evolution is unlikely 
to have a significant role in any under-density observed in 
the APM survey area. 

We are now in a position to examine the number counts 
over the APM survey area. In Fig. [I] we show counts ex- 
tracted for the ~4000 deg 2 field along with the homoge- 
neous and the Northern and Southern 2dFGRS variable <j>* 
models shown in Fig. 121 The form of the cou nts is in good 
agreement with the B (iBusswell et~aT]|2004t) and .Ks-band 



iFrith et al.l2005ab bright number counts measured over the 
APM survey area, with a deficiency of «25 per cent below 
H — 13. In addition, the form of the counts is similar to that 
of the counts extracted from the 2dFGRS Southern declina- 
tion strip and the corresponding variable <f>* model (this is 
also observed in the B and its-band); this perhaps indicates 
that the form of the local galaxy distribution in the ~600 
deg 2 2dFGRS Southern declination strip is similar to that of 
the much larger APM survey area, with an under-density of 
~25 per cent to z =0.1. However, the 2MASS i/-band counts 
over almost the entire survey (|6| >20°, «27000 deg 2 ) are 
also deficient (as are the i^-band counts), with a relatively 
constant deficit of ~15-20 per cent to H = 14 (Fig El right 
hand panels). 

The low |6| >20° counts raise the question as to whether 
systematic effects are significant, or whether these counts 
are due to real features in the local galaxy distribution, as 
suggested by the agreement between the variable <j>* mod- 
els and corresponding counts in Fig. |3] If the latter is true, 
then the size of the local hole would not only be much larger 
than previously suggested but would also represent an even 
more significant departure from the form of clustering at 
large scales expected in a ACDM cosmology. In the follow- 
ing two sections we address a possible source of systematic 
error using new faint //-band photometry - the model nor- 
malisation. Other possible causes for the low counts are also 
discussed in section 5. 
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Figure 5. Here we examine the uncertainty in our photometric 
calibration with 2MASS. The H-band magnitudes determined by 
2MASS and the residual with our photometry are indicated for 
393 point sources below H = 15.1. The large datapoints indicate 
the mean offset and rms dispersion as a function of magnitude. 
The zeropoint used is indicated by the solid line and is accurate 
to ±0.01 magnitudes at lcr confidence. 



Figur e 6. We co mpare the //-band photometry of the LCIRS 
iChen et alj|2002>) with 2MASS using 438 points sources. As in 
Fig- El the large datapoints indicate the mean offset and rms dis- 
persion as a function of magnitude. The mean offset is -0.28±0.01 
magnitudes at lcr confidence. The zeropoint used in the new data 
presented in this work is indicated by the solid line. 



3 NEW FAINT //-BAND DATA 

3.1 Observations &: data reduction 

Our data were taken during a three night observing run 
in September 2004 at the f/3.5 prime focus of the Calar 
Alto 3.5m telescope in the Sierra de Los Filabres in Andalu- 
cia, southern Spain. The J7-2000 infra-red camera contains a 
2048 x 2048 pixel HAWAII- 2 Rockwell detector array, with 
18.5/im pixels, giving a scale of 0.45" /pixel at the prime fo- 
cus. All observations were taken with the //-band filter. Poor 
weather meant that only just over one night's worth of data 
were usable, and even then conditions were not photometric. 

Our primary objective was to image the William Her- 
schel Deep Field as deeply as possible (the results of which 
are presented in a forthcoming paper), but time was avail- 
able at the start of each night to image several 'random' 
fields for 15 minutes each. These were composed of indi- 
vidual 3 second exposures, stacked in batches of 10 before 
readout. A dithering pattern on the sky with a shift up to 
±25" around the nominal centre was adopted. 

Data reduction was complicated by the fact that both 
the dome and twilight sky flat fields appeared to have a 
complicated out-of-focus pattern of the optical train im- 
printed upon them (probably an image of the top end of the 
telescope). This appeared (in reverse) in the science data if 
these frames where used for flat-fielding. We therefore con- 
structed a master flat field by medianing together all the 
science frames from a particular night. This was then used 
to flat-field all the data. Then, individual running medians 
were constructed from batches of 10 or so temporally ad- 



jacent frames, and these were subtracted from each frame 
to produce a flat, background subtracted image. These were 
then aligned and stacked together (with sigma clipping to 
remove hot pixels). 



3.2 Calibration 

Photometric calibration of the //-band images is obtained 
through comparison with the 2MASS point source cata- 
logue. Fig. |H| shows the 2MASS magnitudes compared with 
our data for 393 matched point sources over the Calar Alto 
field and the William Herschel Deep Field. The zeropoint of 
our data is chosen to match that of the 2MASS objects and 
is accurate to ±0.01 magnitudes. The large datapoints and 
errorbars indicate the mean offset and rms dispersion as a 
function of magnitude. When comparing this data to the 
2MASS number counts at bright magnitudes it is important 
to note that the 2MASS point source catalogue includes a 
maximum bias in the photometric zeropoint of <2 per cent 
around the sky (see the 2MASS website). 



3.3 Star/Galaxy separation 

We use the Sextractor software to separate objects be- 
low H — 18; for this magnitude limit, the associ- 
ated STAR_CLASS parameter provides a reliable indi- 
cator of stars and galaxies. We identify 30.0 per cent 
as galaxies (CLASS_STAR<0.1), 58.9 per cent as stars 
(CLASS_STAR>0.9), leaving 11.1 per cent as unclassified. 
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Table 1. The raw number counts per half magnitude are shown for the new //-band data described in section 3 - the CA field (0.27 
deg 2 ) and WHDF (0.06 deg 2 ) in columns 2 and 3. In addition, we show the counts for the LCIRS fields, the HDFS (0.24 deg 2 ) and CDFS 
(0.16 deg 2 ) in columns 4 and 5, applying the zeropoint offset determined with respect to 2MASS in section 4.1. The total number count 
per deg 2 for all fields combined (0.69 deg 2 ) is shown in column 6 along with the homogeneous pure luminosity evolution prediction of 
IMetcalfe etaD <2005l) in column 7. The faintest magnitude bin for the CA field is slightly smaller (0.21 deg 2 ) than at brighter magnitudes; 
the combined solid angle for the faintest bin in column 6 is therefore 0.66 deg 2 . 



4 FAINT //-BAND COUNTS 

4.1 Comparison with the LCIRS 

Before determining number counts for the new //-band data 
described in section 3, we first examine the ph otometry of 
the L as Campanas Infra-Red Survey (LCIRS; IChen et alJ 
2002). The published data covers 847 arcmin 2 in the Hubble 
Deep Field South (HDFS) and 561 arcmin 2 in the Chandra 
Deep Field South (CDFS); the combined solid angle (0.39 
deg 2 ) represents the largest /f-band dataset for 14<//<20. 
The associate d number coun t s are ~15 per cent below the 
homogeneous IMetcalfe et al.l <l2005l) predictions at H = 18 
(see Fig. [TJ. This is significant, as if the model normalisa- 
tion was altered to fit, the deficiency in the 2MASS counts 
at bright magnitudes (Fig. [3J would become much less se- 
vere. However, various other surveys show higher counts, 
although over much smaller solid angles. With the LCIRS 
data in particular therefore, it is vital to ensure that the 
photometric zeropoint is consistent with the 2MASS data 
at bright magnitudes. 

In Fig. El we compare the LCIRS and 2MASS //-band 
photometry for 438 points sources matched over the HDFS 
and CDFS fields. There appears to be a large offset which 
is approximately constant for K >12. Using point sources 
matched at all magnitudes, we determine a mean offset of 
-0.28±0.01 magnitudes; this is robust to changes in the mag- 
nitude range and is consistent over both the HDFS and 
CDFS fields. 

4.2 New //-band counts 

In Fig. |7| we show counts determined for the new //-band 
data described in section 3, the 0.27 deg 2 CA field and the 
0.06 deg 2 WHDF (see also tableQ. Both sets of counts are 
in excellent agreement with the pure lumin osity evolution 
and n o evolution homogeneous predictions of Metcalfe et al. 
(2005). In addition we show LCIRS counts determined in 
the 0.24 deg 2 HDFS and 0.16 deg 2 CDFS, applying the 
0.28 magnitude zeropoint offset determined with respect to 
2MASS in section 4.1. The associated counts are also in ex- 
cellent agreement with the Metcal fe et alJ J2005) models at 
all magnitudes. 

In Fig.|H| we show counts determined from our data and 



the LCIRS combined, with a consistent zeropoint applied as 
in Fig. |7| We estimate the uncertainty arising from cosmic 
variance using field-to-field errors, weighted by the solid an- 
gle of each fie ld. These combined co unts are in good agree- 
ment with the lMetcalfe et al] 1120051) models, particularly at 
fainter magnitudes where the dispersion in the counts aris- 
ing from cosmic variance appears to be small. We perform 
least squares fits between these counts and the pure lumi- 
nosity evolution model; in the magnitude range 14< H <18 
we find a best fi t normalisation o f 1.09 51° n^, where 1.0 cor- 
responds to the lMetcalfe et ail ll2005l) normalisation shown 
in Fig. |H] Varying the fitting range does slightly alter the re- 
sult; in the range 16< H <18 we find a best fit normalisation 
of 1.061+0 033 for example. 



5 DISCUSSION 

In the previous sections, bright //-band number counts from 
2MASS were determined over the APM survey area (~ 4000 
deg 2 ) and almost the entire 66 per cent of the sky (|&| >20°, 
~27000 deg 2 ), along with faint counts to H = 18 over a 
combined solid angle of 0.69 deg 2 applying a zeropoint con- 
sistent with 2MASS. The bright //-band number counts over 
the APM survey area are extremely low (« 25 per cent at 
H — 13) with respect to homogeneous predictions, and re- 
produce the form of the bright c ounts observed in the opti- 
cal B-band (Busswell et al. 2004]) and the near infrared K s - 
band (Frith ct al. 2005a). Previous work has suggested that, 
if due solely to local large-scale structure, these low counts 
would be at odds with the form of clustering expected in a 
ACDM cosmology. In addition, the bright Zf-band |6| >20° 
counts were also found to be low. In the following section, 
various possible causes for these low counts are examined. 



5.1 Model normalisation 

The normalisation of number count models may be deter- 
mined by fixing the predicted to the observed number of 
galaxies at faint magnitudes. The magnitude range at which 
this is done should be bright enough to avoid large uncertain- 
ties in the evolutionary model while faint enough such that 
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Figure 7. Hero wc show //-band galaxy number counts for the 
two separate fields observed in this work, the Calar Alto field (CA 
field; 0.27 deg 2 ) and the William Herschel Deep field (WHDF; 
0.06 deg 2 ). We also show num ber counts determ ined for the two 
separate fields of the LCIRS iChen et alj|2002f) situated in the 
Hubble Deep Field South (HDFS; 0.24 deg 2 ) and Chandra Deep 
Field South (CDFS; 0.16 deg 2 ), subtracting 0.28 magnitudes in 
each case in order to bring the LCIRS and 2MASS zeropoints (and 
hence also the CA field and WHDF zeropoints) into agreement. 
We also show bright number counts extracted from 2MASS for 
the APM survey area and for |6| >20° as shown in Fig. |3 The 
models are indicated as in Fig. In the lower panel, we divide 
through by the pure luminosity evolution homogeneous prediction 
as in Figs. |3 and El At faint magnitudes, wc indicate the Poisson 
uncertainty in each bin. We omit Poisson errors on the bright 
counts for clarity (see Fig.|3]for these). We discuss the uncertainty 
in the counts arising from cosmic variance in section 5. 



large fluctuations in the counts arising from cosmic vari- 
ance are expected to be small. Near infrared wavelengths 
are expected to be insensitive to luminosity evolution at 
bright magnitudes, making the //-band particularly useful 
for such analysis. Of vital importance when determining the 
model normalisation is that when making comparisons be- 
tween faint and bright counts, the zeropoints are consistent; 
an offset of a few tenths of a magnitude between the two, for 
example, would be enough to remove the observed anomaly 
in the bright counts over the APM survey area. 

Applying the 2MASS zeropoint to the fa int //-band 
data presented in this work and the LCIRS data jChen et alJ 
2002), covering a combined solid angle of 0.69 deg 2 , it is clear 
that a discrepancy between the bright and faint counts ex- 
ists; the model normalisation used previously, which indi- 
cates low counts below H = 14 over the APM survey area 
(and for |6| >20°), provides good agreement with the faint 
data. In fact, fixing the model to the faint counts implies a 
slightly higher normalisation. This agreement, as indicated 
by the errorbars in Fig. [g] suggests that the discrepancy 
between the bright and faint counts is not due to cosmic 



Figure 8. Here we show the faint //-band data from the two fields 
presented in this work (CA field and WHDF) and the two field s 
published by the LCIRS (HDFS and CDFS: IChen et al.ll200j) . 
applying a zeropoint to the LCIRS data consistent with the bright 
//-band 2MASS data (and hence the CA field and WHDF also), 
as shown in Fig. [7] The errorbars at faint magnitudes indicate the 
field-to-field error, weighted in order to account for the different 
solid angles of each field. Bright //-band counts extracted from 
2MASS for the APM survey area and for [6| >20° are shown as 
previously. In the lower panel, the counts are divided through by 
the pure luminosity evolution homogeneous prediction as before. 



variance in the faint data. To remove the observed deficit in 
the APM survey area counts below H = 14 by renormalising 
the model, requires a deviation from the faint counts of 7.0a 
using the best fit normalisation of 1.0951q'q34 (determined 
for 14 < H < 18). Similarly, renormalising to the |6| >20° 
counts would require a deviation of 7.2cr from the faint data. 

In addition, the model normalisation may also be scruti- 
nised thro ugh comparison with redshift distributions. Fig.|21 
shows the iMetcalfe et alJ J2005I) pure luminosity evolution 
model compared with //"-band n(z) determined through 
a match between 2MASS and the 2dFGRS Northern and 
Southern declination strips. The model predictions appear 
to be consistent with the observations, with relatively ho- 
mogeneous distributions beyond z = 0.1 (1 per cent and 
8 per cent over-dense in the North and South respectively). 
Lowering the model normalisation to fit the bright 2MASS 
number counts would compromise this agreement and im- 
ply large over-densities beyond z = 0.1 (19 per cent and 
27 per cent in the North and South respectively). 

5.2 Galaxy evolution 

A change in amplitude therefore, cannot easily account for 
the discrepancy in the number counts at bright magnitudes. 
However, could an unexpected change in the slope of the 
number count model contribute? In section 2, we examined 
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the consistency of the number counts at bright magnitudes 
with the underlying redshift distribution, assuming a model 
with insignificant levels of luminosity evolution at low red- 
shift. The predictions derived from the observed n(z) were 
in good agreement with the observed number counts indi- 
cating that luminosity evolution at low redshift is unlikely 
to have a significant impact on the form of the counts at 
bright magnitudes. This is supported by the consistency of 
the pure luminosity evolution model with the observed red- 
shift distributions (Fig. strong low redshift luminosity 
evolution produces a tail in the n(z) which would imply 
large deficiencies at high redshift. 

Could unexpectedly high levels of luminosity evolution 
at higher redshifts affect our interpretation of the bright 
counts? If the slope of the homogeneous prediction were to 
increase significantly above H ~ 14 from the evolutionary 
models considered in this paper, then the model normali- 
sation could effectively be lowered into agreement with the 
bright counts. The problem with this is that the number 
counts beyond H fa 14 are consistent with low levels of lu- 
minosity evolution to extremely faint magnitudes (H ~ 26). 
Models with significantly higher levels of luminosity evolu- 
tion above H » 14 would therefore compromise this agree- 
ment. 

Therefore, it appears that relatively low levels of lumi- 
nosity evolution are consistent with number count obser- 
vations to high redshifts. Also, recent evidence from the 
COMBO-17 survey, examining the evolution of early- type 
galax ies using nearly 5000 objects to z fa 1 jBell et alJ 
2004), suggests that density evolution will also not con- 
tribute; <f>* appears to decrease with redshift indicating that 
the number of objects on the red sequence increases with 
time, and so acts contrary to the low counts observed at 
bright magnitudes. This p icture is supported by the K20 
survey dCimatti et al.l2002t) , which includes redshifts for 480 
galaxies to a mean depth of z fa 0.7 and a magnitude limit 
of K a = 20 with high completeness. The resulting redshift 
distribution is cons istent with low leve ls of luminosity and 
density evolution jMetcalfe et alJl2005ft . 

In summary, significant levels of evolution are not ex- 
pected in passive or star forming pure luminosity evolu- 
tion models, although could occur through dynamical evo- 
lution. However, the pure luinosity evolution models of 
iMetcalfe etaH J2005h fit the observed H < 14 n(z) at 
2 > 0.1; it is at lower redshifts that there are fluctuations. In 
addition, these models continue to fit the observed n(z) at 
very high redshift and the number counts to extremely faint 
magnitudes (K fa 23), suggesting that there is little need for 
evolution at z fa 1, far less z< 0.1. Some combination of dy- 
namical and luminosity evolution might be able to account 
for these observations; however it would require fine-tuning 
in order to fit both the steep counts at bright magnitudes 
and the unevolved n{z) at low and high redshifts. 

5.3 Photometry issues & completeness 

The number counts shown in Figs. Q and |H] show bright and 
faint counts with a consistent zeropoint applied. Photom- 
etry comparisons have been made using several hundred 
point sources matched at bright magnitudes. In order to 
check that the applied zeropoints are consistent with the 
galaxy samples, we also compare the 2MASS photometry 



with 24 matched galaxies in the CA field and WHDF and 
16 in the LCIRS samples; we find that the mean offsets are 
— 0.01±0.04 and — 0.32±0.06, consistent with the zeropoints 
determined via the 2MASS point sources. The comparisons 
with the 2MASS point source catalogue (Figs. |K|and|SJl also 
indicate that there is no evidence of scale error in either of 
the faint samples to H fa 16. 

Could the discrepancy between the bright and faint 
counts arise from an under-estimation of the total flux of the 
galaxies? Recall that we make no correction to total mag- 
nitude for the faint data presented in this work; however, 
under-estimating the total flux in the faint data would only 
increase the observed deficit in the counts at bright magni- 
tudes, if the model normalisation is adjusted to fit the faint 
counts. The good agreement between the point source and 
galaxy zeropoints suggests that the estimate for the total 
galaxy flux is comparable in the bright and faint data. At 
bright magnitudes, the 2MASS extrapolated ff-band mag- 
nitudes are used. In the A's-band, this magnitude estimator 
has been shown to be an excellent estimate of the total flux, 
through c omparisons w i th th e total _R" 3 -band magnitude es- 
ti mator of I Jones et alJ 11200411 and the if-band photometry 
of lLovedavH2000l) . 

Another possible contribution to the low counts could 
be high levels of incompleteness in the 2MASS survey. 
As with the possible systematic effects described previ- 
ously, it is differing levels of completeness in the faint and 
bright data which would be important. The 2MASS litera- 
ture quotes the extended source catalogue completeness as 
> 90 per c ent (see the 2MA SS website for example). Inde- 
pendently, iBell et all (I2003T) suggest that the level of com- 
pleteness is high (~99 per cent), determined via comparisons 
with the SDSS Early Data Release spectroscopic data and 
the 2dFGRS. The faint data presented in this work and the 
LCIRS data are likely to suffer less from incompleteness, as 
we cut well below the magnitude limit, are subject to lower 
levels of stellar confusion and suffer less from low resolution 
effects. Incompleteness in 2MASS will therefore affect the 
observed deficit in the bright counts at the < 10 per cent 
level, although the effect is likely to be at the low end of 
this constraint due to incompleteness in the faint catalogues 
and suggestions that the 2MASS extended source catalogue 
is fairly complete. 

5.4 Large-scale structure 

It appears therefore, that the observed deficiency in the 
bright counts might be significantly affected by incomplete- 
ness in the 2MASS extended source catalogue. However, the 
level to which other systematic effects such as the model 
normalisation, luminosity evolution and photometry issues 
appears to be small. The question then is — accounting for 
these various sources of error or uncertainty, are the defi- 
ciencies in the bright i/-band counts over the APM survey 
area and for |6| > 20° still at odds with the expected fluc- 
tuations in the counts arising from local large-scale struc- 
ture in a ACDM cosmology, as suggeste d in previous work 
jBusswell et alJl2004l : iFrith et alJl2005al) ? 

We determine the expected fluctuations in the bright 
number counts due to cosmic variance via ACDM mock 
2MASS catalogue s; these are described in detail in 
IFrith et alJ J2005aft . To recap, we apply the 2MASS selec- 
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tion function to 27 virtually independent volumes of r = 
500 h' 1 Mpc formed from the 3000 3 /i" 3 Mpc 3 ACDM Hubble 
Volume simulation. This simulation has input parameters of 
tt m = 0.3, Q b = 0.04, h = 0.7 and a s = 0.9 ilenkins et alJ 



1998). The mean number density of the counts at the mag- 
nitude limit is set to that of the observed 2MASS density. 

We are now in a position to estimate the significance 
of the observed bright H-band counts. We use the la fluc- 
tuation in the counts expected in a ACDM cosmology (de- 
termined using the 2MASS mocks described above), which 
for the APM survey area is 7.63 per cent (for H < 13) 
and 4.79 per cent (for H < 14), and for |6| > 20° is 
3.25 per cent (for H < 13) and 1.90 per cent (for H < 14). 
In addition we also take into account the uncertainty in the 
model normalisation; we use the best fit normalisation of the 
iMetcalfe et alJ <|2005T) p ure luminosity evoluti on model (a 
factor of 1.095 above the IMetcalfe et all i2005l) model) and 
add the uncertainty of ±3.1 per cent derived from the faint 
//-band counts (presented in Fig. |HJ in quadrature. Regard- 
ing the possible effect arising from survey incompleteness, 
we first assume that the level of incompleteness is compara- 
ble in the faint and bright data; the resulting significance for 
the APM survey area and \b\ > 20° bright counts are shown 
in column 3 of tabled This represents an upper limit on the 
significance since we have effectively assumed that there is 
no difference in the incompleteness between the bright and 
faint datasets. In column 4 of table |5] we assume that there 
is a difference in the completeness levels in the faint and 
bright data of 10 per cent. This represents a lower limit on 
the significance (assuming that there are no further signifi- 
cant systematic effects) , since we assume that the complete- 
ness of the 2MASS extended source catalogue is 90 per cent 
(the lower limit) and that there is no incompleteness in the 
faint data. 

Therefore, assuming a ACDM cosmology, it appears 
that the observed counts over the APM survey area might 
be in line with a rare fluctuation in the local galaxy dis- 
tribution. However, the counts over 66 per cent of the sky 
(\b\ > 20°) suggest a deficiency in the counts that are at 
odds with ACDM, even accounting for a 10 per cent incom- 
pleteness effect and the measured uncertainty in the best fit 
model normalisation. 



6 CONCLUSIONS 

We have presented new //-band photometry over two fields 
with a combined solid angle of 0.30 deg 2 to H «19. The 
zeropoint is chosen to match that of the 2MASS photometry 
at the bright end and is accurate to ±0.01 magnitudes. In 
addition we have examin ed the faint //-band data of the 
LCIRS llChen et al.ll2002l) which covers two fields with a 
combined solid angle of 0.39 deg 2 to H «20. The zeropoint of 
this data appears to be offset from the 2MASS photometry 
by 0.28±0.01 magnitudes. Applying a consistent zeropoint, 
the faint counts determined from the new data presented in 
this work and the LCIRS are in go od agreement with th e 
pure luminosity evolution model of IMetcalfe et alJ p005), 
with a best fit normalisation a factor of 1.095lg 034 higher. 

In contrast, the bright //-band counts extracted from 
2MASS over the «4000 deg 2 APM survey area around the 
southern galactic pole are low with respect to this model, 



Field 




Significance 
(no incompleteness 
correction) 


Significance 
(assuming 10 per cent 
incompleteness) 


APM 


13.0 


3.7a 


2.5a 


APM 


14.0 


4.2a 


2Aa 


|6| > 20° 


13.0 


6.1a 


3.8a 


6 > 20° 


14.0 


6.8a 


4. Oct 



Table 2. Here we show the significance of the //-band 2MASS 
counts extracted for the ss4000 deg 2 APM survey area and for 
\b\ > 20°, for H < 13 and H < 14. In each case we determine the 
expected cosmic variance using a ACDM form to the large-scale 
power determined via mocks constructed from the Hubble Vol- 
ume simula tion. In addit ion we use the best fit normalisation of 
the IMetcalfe et alJ 120051) pure luminosity evolution model deter- 
mined at faint magnitudes (a factor of 1.095 above the original 
normalisation) and add the uncertainty on this (±3.1 per cent) in 
quadrature to the expected cosmic variance. In the third column 
we use the observed counts as shown in Figs . HI 171 and IS1 in the 
fourth column we account for an upper limit on the incomplete- 
ness in the 2MASS extended source catalogue of 10 per cent; the 
level to which this will affect the significance is likely to be lower 
due to incompleteness in the faint data. 



corroborating previous results over th is area in the optical 
/?-band and near infrared ZG-band feusswell et al.l 12004 
iFrith et alJl2005al) . In addition, the counts extracted for al- 
most the entire survey, covering 66 per cent of ths sky, are 
also low with a deficit of 15 — 20 per cent to H = 14. Im- 
portantly, this discrepancy does not appear to be due to ze- 
ropoint differences between the faint and bright data or un- 
certainty in the model normalisation set by the faint counts. 

We have investigated various possible sources of sys- 
tematic error which might affect this result: The counts are 
consistent with low levels of luminosity and density evolu- 
ti on, as predicted by t he pure luminosity evolution model 
of lMetcalfe et all <l2005l) , to extremely faint magnitudes (see 
Fig. 0. Also, the photometry appears to be consistent be- 
tween the faint and bright galaxy data using a zeropoint 
applied via comparisons between point sources. However, 
differing incompleteness in the bright and faint galaxy sam- 
ples might have a significant impact; completeness in the 
2MASS extended source catalogue is > 90 per cent. 

Finally, we determined the expected cosmic variance in 
the bright number counts from ACDM mock 2MASS cat- 
alogues. Allowing for the model normalisation uncertainty 
determined from the faint counts, and using an upper limit 
on the incompleteness in the 2MASS galaxy sample, the 
deficiency in the counts over the APM survey area repre- 
sent a rare («1 in 100) fluctuation in a ACDM cosmology. 
However, the low //-band counts for |6| > 20° suggest that 
this deficiency might extend over the entire local galaxy dis- 
tribution; allowing for incompleteness and the model nor- 
malisation uncertainty as before, this would represent a 4cr 
fluctuation (<1 in 10 000) in the local galaxy distribution, 
and would therefore be at odds with the expected form of 
clustering expected in a ACDM cosmology on large scales. 
The increase in faint near infrared data from the UK In- 
frared Deep Sky Survey (UKIDSS) should help to resolve 
this issue. 
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